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FOREWORD 


This document describes a method for transforming coordinates, developed 
for use in studies of high latitude energetics. The work was conducted by the 
Electro-Mechanical Division, Northrop Corporation, Huntsville, Alabama, for 
the National Aeronautics and Space Administration, George C. Marshall Space 
Flight Center, Aero-Astrodynamics Laboratory, under contract NAS8-20082, 

Appendix A-l, schedule order 32. This program was under the direction of the 
Space Environment Branch, with Mr. R. E. Smith as NASA/MSFC Technical Coordinator. 
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ABSTRACT 


Based upon a mathematical model of contours of constant corrected geomagnetic 
latitude in a polar projection of geographic coordinates, analytic equations 
are developed for converting geographic coordinates to corrected geomagnetic 
coordinates and vice versa* The equations have been programmed for use on 
a small computer. Values generated by these programs are compared with 
published values; accuracies should be sufficient for many studies of high 
latitude geophysical phenomena. This treatment is restricted to the northern 
hemisphere. 
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Section I 

INTRODUCTION 

The purpose of this memorandum is to document the derivation and programm- 
ing of a set of empirical analytic transformations between geographic coordinates 
and corrected geomagnetic coordinates. The need for such transformations arises 
in investigations of the interactions between high latitude geophysical pheno- 
mena (e.g., aurorae, polar magnetic substorms, etc.) and the neutral atmosphere. 
Because high latitude phenomena are controlled by the geomagnetic field, studies 
of these phenomena employ geomagnetic coordinates. For the neutral atmosphere, 
on the other hand, dynamic effects associated with the rotation of the earth 
are more important; hence, geographic coordinates based on the earth 1 s axis of 
rotation, are the natural frame of reference for such investigations. In order 
to employ the results of studies in both of these areas together on a large 
scale, transformations between the coordinate systems are required. 

Over the years, several systems of geomagnetic coordinates have been 
employed, those based on the centered and eccentric dipole approximations of 
the geomagnetic field being most prominent. It has recently become apparent 
that the centered dipole approximation is inadequate for quantitative studies. 
Hultqvist (ref. 1) made a detailed study of the geomagnetic field lines using 
the first five terms in the spherical harmonic expansion of the field. Based 
on these results Hakura (ref. 2) has developed a system of geomagnetic 
coordinates and prepared tables and maps relating them to geographic coordinates. 
He has designated these ’’corrected geomagnetic coordinates”; this terminology 
is followed here. A number of investigators have employed corrected geomagnetic 
coordinates in detailed morphological studies of auroral and geomagnetic activity 
(e.g., refs. 3, 4, 5). Consequently, this system is considered to be the most 
useful for relating such phenomena to effects in the neutral atmosphere. 

Hakura f s (ref. 2) corrected geomagnetic coordinates are obtained from 
tabulated data prepared by Hultqvist (ref. 1); and tabulated data are the 
result of his computations. Direct use of Hakura *s coordinates on a computer 
would require extensive tables relating corrected geomagnetic coordinates to 
geographic coordinates and an interpolation routine to provide the necessary 
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coordinate transformations. As an alternative to such a procedure, a set of 
empirical analytic equations is developed to provide a relatively fast and 
simple means of accomplishing these transformations at a small loss in 
accuracy. This is accomplished by formulating a mathematical model of contours 
of constant corrected geomagnetic latitude in the geographic coordinate system. 
Based on this model, analytic equations are obtained relating the two coordinate 
systems. Details of the model and the derivation of these equations are con- 
tained in Section II. In Section III, computer programs are described and 
errors arising from the transformations are discussed. Section IV contains a 
summation. Computer programs and tables relating the two coordinate systems 
are presented in the appendixes. 
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Section II 

TRANSFORMATION EQUATIONS 


2.1. METHOD 

The approach taken is intuitive, suggested by the resemblance of the 
contours of constant corrected geomagnetic latitude, when drawn on an equidistant 
polar projection of the geographic coordinate (G C) system, to a family of ellipses 
(see Figure 2-1) . Since these contours are circles in the corrected geomagnetic 
coordinate (CGC) system, one step in the conversion of CGC to GC is the trans- 
formation of circles to ellipses. The remaining operations are the usual 
translation of the origin and rotation of the coordinate frame to the established 
reference direction. 

In the conversion of GC to CGC, it has not been feasible to invert the 
above transformation equations. Rather a similar procedure has been followed 
in reverse sequence, making use of the same parameters and some of the same 
equations. 


All calculations are carried out in polar coordinates, r, <j> , in which r 
represents colatitude. Where constant parameters are concerned, only the 
northern hemisphere is considered. 


2.2 CORRECTED GEOMAGNETIC COORDINATES TO GEOGRAPHIC COORDINATES 

As previously indicated, circles of constant latitude in the CGC system 
transform to ellipses in the GC system. Figure 2-2 shows the family of 
ellipses which are taken to approximate the CGC latitude contours in Figure 2-1. 
These ellipses are scaled (same eccentricity) from the 70 degree ellipse. 


Let r, <J> represent coordinates of the polar coordinate system shown in 
Figure 2-2. The equation of an ellipse in polar coordinates is 

2 


(r T V 


2 k 2 

a b 


2 . 2 
a sm 


2 2 
+ b cos 


( 1 ) 


where a and b are the semi-axes of the ellipse. 
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Figure 2-1. CONTOURS OF CONSTANT CORRECTED GEOMAGNETIC LATITUDE IN A NORTH POLAR 
PROJECTION OF THE GEOGRAPHIC COORDINATE SYSTEM. (PREPARED FROM TABLES 
IN REFERENCE 2) 
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Figure 2-2. MODEL CONTOURS OF CONSTANT CORRECTED GEOMAGNETIC LATITUDE IN A POLAR 
PROJECTION, CENTERED ON THE CORRECTED GEOMAGNETIC NORTH POLE AND 
DIRECTION REFERENCED TO THE MAJOR AXIS 





M-765 


NORTHROP 

Huntsville 


Consider a circle of constant latitude, represented by colatitude r. The 
equation of this circle is 

r = constant (2) 

Equation (2) can be put into the form of equation (1) by writing 


2 

r = 


2 u 2 
a b 


2 . 2 ■ u 2 2 
a sin a 4* b cos a 


(3) 


where a is the (constant) value of <J> for which the circle (r) and the ellipse 
(r f ) intersect (see Figure 2-2). Taking the ratio of equations (1) and (3) 
gives 


Cl !) 2 = 
r 2 


sin a + b cos a 

2 2 2 

a sin <f> + b cos <|> 


(4) 


The right side of equation (4) can be rewritten in terms of the eccentricity, 
which is defined as 

e = \j l-b 2 /a 2 

With this parameter it is easy to put equation (4) into the form 

f 2 . 2 2 

> r v _ 1 - e cos a 

~ 2 2 

1 - e cos <J> 

Hence circles transform to ellipses by the equation 

~ 2 2 
1 - e cos a 


i 2 2 

1 - £ COS 


1/2 


( 5 ) 


Now a corresponding coordinate <p 9 must also be determined (it was found 
empirically that merely setting <J> f = <f> gave unacceptable results). A heuristic 
means of picturing the alteration of the <f> coordinate (and one which leads to 
acceptable results) is to imagine the circle of constant latitude being distorted 
into an ellipse by displacements (of points on the perimeter) which are normal 
to the circle and outward at the major axis and inward at the minor axis. For 
each quadrant a constant (small) displacement vector is established by the 
vector sum of the components at the adjacement semi-axes. The <f> component of 
this displacement vector then gives the required modification for each value of 
<p in that quadrant. 
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Explicitly, let the polar coordinates r, <j> specify the vector r, i.e. 

r = r(r,<j>) (6) 

Then r' is represented as 

r' = r' (r + Ar, <t> + A<J>) (7) 

The displacement vector is taken to be 

Ar = r - r (8) 

-y ■> 

where Ar is treated as having small magnitude. Now Ar may be resolved along 
the major and minor axes as 


Ar = i Aa cos <j> + j Ab sin <f> (9) 

where i and j are unit vectors along the major and minor axes respectively and 
Aa and Ab are the corresponding (empirically known) displacements along these 
directions (opposite signs on opposite sides of the circle). In terms of 
polar coordinates, Ar may be written 

Ar = Ar r + r 1 A<j> <j> (10) 

where r and <j> are the unit vectors in the r and < p directions, with Ar and r T A$ 
the corresponding displacements. It will be noted that r ! is used instead of 
r in equation (10). If Ar were sufficiently small this distinction would not 
matter. In fact, Ar is not that small and r 1 gives slightly better results. 
The unit vectors in equations (9) and (10) are related by 


r = i cos (j) + j sin (p 

^ A /V 

<j> = -i sin <j> + j cos 

i 


(ID 
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These may be inverted to give 


i = r cos (p - <f> sin <f) 

/v A A 

j = r sin <p + <f> cos <f> 


( 12 ) 


Placing these results into equation (9) and substituting that result into equation 
(10) gives 


Aa cos^cf) r - ^ sin 2<p <f> + Ab sin^<j> r + sin 2<J> <p 


= Ar r + r T Ac#) 


(13) 


Since only A<J> is required, only <j> components need be considered; hence 


A* = 

r T 


Ab 

2 


Aa 

2 


sin 2<p 


(14) 


Because of the scaling which is used in generating the ellipses, both Aa and Ab 
are proportional to r; therefore 


Ab - Aa = cr 


With this substitution equation (14) becomes 


A<f> = f ( fr ) sin 2<t> 


Taking the ratio from equation (5) gives 


A(j> = - .13 


i 2 2 .-,1/2 

1 - e cos 


i 2 2 

1 - £ cos a 


sin 2<f> 


(15) 


where - .13 has been empirically determined. The coordinate cp' corresponding 
to r' is therefore given by 


= * 


+ A(f> 

- .13 


i 2 2 , 

1 - e cos <p 

. 2 2 

1 - e cos a 


1/2 

sin 2(f) 


(16) 
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corresponds to <J> in the above treatment. 

The origin must now be shifted from the corrected geomagnetic north pole 
to the north geographic pole and a new reference direction corresponding to the 
GC system employed. The relevant geometry is shown in Figure 2-3. From the 
standard relationships in a plane triangle 

(r") 2 = R 2 + r' 2 - 2Rr' cos r' (22) 

where 

T' E 0' - 3 (23) 

and R is the latitudinal distance between the poles in degrees geographic 
latitude. Another relation gives 

rj = _R (24) 

sin y sin 4 > 
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Figure 2-3. GEOMETRIC RELATIONSHIPS IN CHANGING POINT AND DIRECTION OF REFERENCE 
FROM THOSE OF THE CORRECTED GEOMAGNETIC COORDINATE SYSTEM TO THOSE 
OF THE GEOGRAPHIC COORDINATE SYSTEM 
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where the angles are defined by Figure 2-3. From the geometry in Figure 2-3, 
it is clear that 


<P = TT - (r ' + y) 

0" = tt + e - (a + y) 


(25) 

(26) 


where 6 and A are known constant angles. Substituting equation (25) into (24) 
gives 


r 

R 


sin y 


sin [tt - (r' + y)] 


(27) 


With an expansion of the denominator on the right and rearrangement of terms, 
equation (27) may be put into the form 


tan y = 


sin f* 


R/r' - cos T' 


(28) 


so that 


Y = ARCTAN 


sin f' 


R/r' - cos T' 
Substituting this result into equation (26) gives 

6" = tt + 6 - A - ARCTAN 


(29) 


sin T' 


(30) 


LR/r' - cos r' J 

In summary, geographic latitude 0 and longitude 0 are given in terms of 
corrected geomagnetic coordinates Q^, 0^ by the following equations: 


)" = 


CD 

II 

VO 

O 

0 

1 

r" 



(31) 

$ = e" 




(32) 

[R 2 + r' 2 

2r* R cos 

r* ] 1/2 

(22) 

TT + 6 “ 

A - ARCTAN 

- 

sin r' 

(30) 


R/r' - cos T' J 
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. i _ 


1 - 


= r 


2 2 -| 1/2 

e cos a 


L 1 - 


2 2 
e cos T 


r ' = e * - 3 


(19) 

(23) 


r 

0' 


e 

e 


6 



2 2 ' 
e cos r 


2 2 
e cos a 


1/2 

sin 2T 


( 21 ) 

( 20 ) 


r = 90° - 0 

c 


(17) 


0 = $ (18) 
c 

where a, 8, A, e, and R are empirically determined constants whose values are 
listed in Table 2-1. 


Table 2-1. VALUES OF TRANSFORMATION PARAMETERS 


PARAMETER 

VALUE 

a 

55° 

8 

o 

o 

A 

o 

O 

e 

.5913 

R 

9.5° 


2.3 GEOGRAPHIC COORDINATES TO CORRECTED GEOMAGNETIC COORDINATES 

Although it has not been feasible to invert the preceding equations, the 
same parameters are used in a similar manner to obtain the inverse trans- 
formation. The pertinent geometry is given in Figure 2-4, where geographic 
latitude and longitude 0, $ are related to the variables shown by equations 
(31) and (32). 


Translating the origin from the north geographic pole to the corrected 
geomagnetic north pole gives 
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Figure 2-4. GEOMETRIC RELATIONSHIPS IN CHANGING POINT AND DIRECTION OF REFERENCE 
FROM THOSE OF THE GEOGRAPHIC COORDINATE SYSTEM TO THOSE OF THE 
CORRECTED GEOMAGNETIC COORDINATE SYSTEM 
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(r') 2 = (r") 2 + R 2 - 2r" R cos y 


(33) 


From Figure 2-4 


y = 2 tt - 0" - A - (tt - 6) 
= it - 0" + 6 - A 


(34) 


With equation (34) , equation (33) may be rewritten 


r ' = [(r") 2 + R 2 + 2r" R cos (0" - 6 + A)] 1 ^ 2 (35) 


Also from Figure 2-4 and the trigonometry of plane triangles 


where 


sin T sin ip 


ip = n - (T + y) 


(36) 

(37) 


Substituting equation (34) into (37) gives 


^=e"-e + A- r 


(38) 


hence 


r" = sin T 

R “ sin (0" - 6 + A - T) 


The denominator on the right may be expanded and the equation solved for F to 
obtain 


f = ARCTAN 


sin (0" - 6 + A) 

R/r" + cos (0" - 6 + A) 


(39) 


The reference direction of the coordinate system must now be aligned with the 
zero meridian of the CGC system. This is accomplished through (see Figure 2-4) 

0' = r + 6 (40) 
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Finally the coordinate system must be distorted so that ellipses become 
circles. The results of the previous subsection may be used. Combining 
equations (35) and (19) gives 

\ r 2 2 “1 |l/2 

r = [ (r") 2 + R 2 + 2r" R cos (6" - 3 + A) ] i 1 ~ £ 2 COS 2 ~ r 

1 1 “ e cos aj i 


(41) 


Similarly, combining equations (40) and (20) gives 

, 1/2 


0 = r + 6 + .13 


2 2 \ 

1 - £ cos T ! sin 2T 

i 2 2 

1 - e cos a . 


(42) 


Summarizing, the conversion of geographic to corrected geomagnetic coordinates 
may be effected per the following equations: 


0 - 90° - r 

c 

$ = e 

c 

r = [ (r") 2 + R 2 + 2r" R cos (0" - 3 + A)] 1 1 


(17) 


(18) 


2 2 -j ] 1/2 

- e cos f ! ; 


2 2 
e cos a 


(41) 


1 - 


e = r + 3 + .13 


, 2 2 

1 - £ COS r 

1 2 2 

1 - e cos a 


1/2 

sin 2r 


T = ARCTAN 


sin (9" - 3 + A) 


R/r" + cos (6" - 3 + A) 


(42) 


(39) 


r" = 90° - 0 


0 " = $ 


(31) 

(32) 


where R, a, 8» A, e are known empirical constants whose values are listed in 
Table 2-1. 
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Section III 
RESULTS 


3.1 COMPUTER PROGRAMS 

Equations for transforming corrected geomagnetic coordinates to geographic 
coordinates have been programmed in Fortran for use on the IBM 1130 computer. 

This program, designated MC2GC, together with computed tables is presented in 
Appendix A. A similar program, designated GC2MC, has been prepared for 
transforming geographic to corrected geomagnetic coordinates. Appendix B 
contains this program and computed tables. Computer time required for each 
program to produce the tables shown, including compilation time, was less 
than 6 minutes, determined primarily by the speed of the printer. 

Both programs have been developed with an input option which permits input 
coordinates to be either read in or generated internally. The latter is 
useful for synoptic studies, while the former is useful for fixed or instaneous 
points of interest, such as geophysical stations, satellite passes, or boundaries 
of the auroral oval. Tables presented in the appendixes have been produced 
through the synoptic option, covering latitudes 50 degrees at 5-degree 
intervals and all longitudes at 10-degree intervals. 

Minor program modification will permit changes of coverage and intervals. 

In addition to a switch variable designating the input option, required input 
includes the transformation parameters listed in Table 2-1, and, unless the 
synoptic option is designated , the input coordinates together with the number 
of coordinates pairs for a given run. 

3.2 ERRORS 

It is useful to establish a common unit in which all errors are to be 
expressed, so that deviations in both latitude and longitude for both trans- 
formations have the same relative significance. For these purposes, it is 
convenient to take a degree of geographic latitude as such a unit. This unit 
is proportional to distance on the surface of the (reference sphere) earth; 
that is, along a geographic meridian, arc length is given by 

S = y AG I (43) 

$ = constant 
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where 


it R 

360° 6 


(44) 


with the radius of the earth and A0 given in degrees. In units of y then, 
distance on the surface of the earth (arc length) is given by 


S = A0 | 

0 = constant 


(45) 


Similarly it is easy to see that arc length along a circle of constant geographi 
latitude is given by 


S = y cos© A$ | 

0 = constant 


So in units of y, 

S = cos0 A$ | = (A$) f (46) 

0 = constant 

When A$ and A0 are expressed in degrees, both (A$) f and AO represent distances 
in the same way; they may therefore be compared with one another. 


The CGC system has an additional consideration. Since translations and 
rotations preserve distances, they introduce no alterations into equations 
(45) and (46). However, this is not true for the transformations of circles 
to ellipses and vice versa. In the CGC system arc length along a meridian of 
longitude is still proportional to a change in latitude, as in equation (45), 
but the constant of proportionality varies with longitude. This variation is 
just that expressed by equation (19). Hence, from equations (17), (18), (19), 
(21) , and (45) 


S 


A0 


1 - cos^ ($> - 8) 

c 


1 - 


2 2 

e cos a 


(A0 c ) 1 (47) 




constant 
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where jj is taken as the unit of length. From equation (46) 

S = cos 0 A<h I l (A^ ) f (48) 

c c 1 c 

0 = constant 

c 

where it is assumed that the variation of A0 in equation (20) over the interval 
of A$ c is negligible. The deviations A0, (A$) f , (A0^) 1 , and (A$^) f , as defined 

by equations (45) - (48) , are all proportional to distance on the earth in the 
same way, and on this basis may be compared with one another. 

The corrected geomagnetic coordinates of Hakura (ref. 2) are taken as the 
standard for comparison. Errors are defined as deviations of present results 
from those in that reference. The extensive tables presented there are used 
as reference values with coordinate conversions in both directions being 
checked against them. 


For the transformation of geographic to corrected geomagnetic coordinates, 
comparison between the tables in Appendix B and the tables in reference 2 is 
direct and convenient. Maximum positive and negative deviations are computed 
according to 


AX 5 Max 


(X 


computed 


X Hakura^ 


(49) 


Table 3-1 contains maximum positive and negative deviations in corrected 
geomagnetic latitude (A0^) T for each given geographic latitude shown. Also 
shown is the longitude at which the maximum deviation occurs and the standard 
values of the corrected geomagnetic coordinates at that point. Similar 
information for (A#^) T , obtained in the same manner, is presented in Table 
3-2. 


From Tables 3-1 and 3-2 some general conslusions can be drawn. Errors 
appear to increase with decreasing latitude, so that the latitude limits for 
which the transformations are useful depend upon the inaccuracies which can be 
tolerated. Corrected geomagnetic latitude and longitude errors are comparable 
for given geographic latitude, although the maximum deviations for both do not 
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Table 3-1. MAXIMUM DEVIATIONS OF COMPUTED VALUES OF CORRECTED 

GEOMAGNETIC LATITUDE FROM STANDARD VALUES ALONG CURVES 
OF CONSTANT GEOGRAPHIC LATITUDE 


0 

(A0 C )'* 

$ 

0 ** 
c 

$ ** 
c 

50° 

+ 

2.1° 

140° 

43.9° 

210.0° 



.3° 

330° 

54.9° 

56.0° 

55° 

+ 

1.7° 

150° 

48.7° 

217.8° 


- 

.4° 

330° 

59.8° 

58.3° 

O 

O 

LTD 

+ 

1.5° 

150° 

53.8° 

216.7° 


- 

.2° 

330° 

64.4° 

60.9° 

O 

LO 

LO 

+ 

1.1° 

160° 

59.1° 

222.3° 



.4° 

290° 

76.6° 

11.0° 

O 

O 

+ 

.9° 

160° 

64.2° 

219.3° 


- 

.4° 

270° 

80.9° 

327.6° 

75° 

+ 

.7° 

170° 

69.8° 

221 .0° 


- 

.3° 

280° 

85.7° 

353.1° 

0 

o 

00 

+ 

.5° 

170° 

74.6° 

212.7° 


- 

.1° 

60° 

73.9° 

144.8° 

00 

cn 

o 

+ 

.4° 

190° 

79.9° 

202.6° 


* Defined by equations (47) and (49) 

** Standard values - from Hakura (ref. 2) 
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Table 3-2. MAXIMUM DEVIATIONS OF COMPUTED VALUES OF CORRECTED 
GEOMAGNETIC LONGITUDE FROM STANDARD VALUES ALONG 
CURVES OF CONSTANT GEOGRAPHIC LATITUDE 


0 

(A<t> c )'* 

$ 

0 ** 
c 

■ 0 ** 
c 

50° 

+ 1.4° 

330° 

54.9° 

56.0° 


- 2.2° 

230° 

53.9° 

285.3° 

O 

LO 

LO 

+ .9° 

330° 

59.8° 

58.3° 


- 1.5° 

230° 

59.0° 

282.9° 

60° 

+ .7° 

330° 

64.4° 

60.9° 


- 1.0° 

230° 

63.9° 

279.8° 

65° 

+ .4° 

320° 

71.2° 

54.0° 


- 1.1° 

130° 

59.3° 

198.6° 

O 

0 

+ .3° 

200° 

67.9° 

246.4° 



130° 

64.1° 

196.8° 

75° 

+ .5° 

290° 

84.9° 

26.8° 


- .7° 

130° 

68.8° 

194.1° 

00 

o 

o 

+ .5° 

230° 

81 .3° 

239.9° 


- .4° 

130° 

73.3° 

190.0° 

0 

LO 

00 

+ .2° 

240° 

83.7° 

203.3° 


- .1° 

130° 

77.7° 

183.0° 


* Defined by equations (48) and (49) 

** Standard values - from Hakura (ref. 2) 
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necessarily occur in the same direction (longitude). Not indicated in the 
tables but apparent in the comparison is the fact that within 1-2 degrees 
latitude of the geomagnetic pole (CGC system) longitude errors become large in 
terms of degrees, although not in terms of distance as expressed by (A30 T due 
to the cos 0^ factor. In essence, the longitude conversion breaks down near the 
pole; but since the latitude conversion holds well in this region, problems of 
physical location are not very serious. 


For the conversion of corrected geomagnetic coordinates to geographic 
coordinates, the tables of Hakura (ref. 2) are not convenient for comparison. 
Nevertheless, they can be used because they do establish a correspondence 
between coordinates in the two systems. In order to check the accuracy of this 
transformation, corrected geomagnetic coordinates of points along contours of 
constant geographic latitude were taken from reference 2 and read into program 
MC2GC. The computed geographic coordinates were then compared with those from 
the reference, to which the input corrected geomagnetic coordinates corre- 
sponded. Tables 3-3 and 3-4, for errors in geographic latitude and longitude 
respectively, were prepared in this way, according to the notation of equation 
(49). The correct values of all coordinates at the point where the maximum 
deviation occurs are also included. It should be noted that to be strictly 
analogous with Tables 3-1 and 3-2, contours of constant corrected geomagnetic 
latitude should have been used in determining the maximum deviations. Since 
this was impractical, given the tables of reference 2, trends in Tables 3-3 and 
3-4 are not apparent, unlike Tables 3-1 and 3-2. *For present purposes, it is 
sufficient to note that errors for this transformation are comparable in 
magnitude to those of the previous one. 


To supplement Table 3-3, Figure 3-1 has been prepared. This figure compares 
computed contours of constant corrected geomagnetic latitude (solid line) with 
those of Hakura (ref. 2) (dashed line) in the geographic coordinate system. 

Since geomagnetic latitude is generally the coordinate of importance when the 
geomagnetic field is important, corrected geomagnetic longitude (computed) is 
indicated only at the intersection with the (computed) latitude contours. Figure 
3-1 clearly shows the "distribution of error 11 in corrected geomagnetic latitude 
arising from the transformation. From a comparison of Tables 3-3 and 3-4 it 
can be seen that longitudinal errors are not necessarily distributed in the same 
manner, particularly at higher latitudes. 
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Table 3-3. MAXIMUM DEVIATIONS OF COMPUTED VALUES OF GEOGRAPHIC 
LATITUDE FROM STANDARD VALUES ALONG CURVES OF 
CONSTANT GEOGRAPHIC LATITUDE * 


0 

A0 

$ 

0 ** 
c 

$ ** 
c . 

40.0° 

+ .5° 

310.0° 

50.7° 

30.9° 


- 2.6° 

265.0° . 

51.4° 

328.2° 

50.0° 

+ .8° 

315.0° 

58.8° 

40.3° 


- .9° 

215.0° 

50.5° 

271.5° 

60.0° 

+ 1.1° 

45.0° 

55.5° 

119.5° 


- .8° 

175.0° 

54.5° 

235.6° 

70.0° 

+ .9° 

50.0° 

65.1° 

128.2° 


- .4° 

190.0° 

66.5° 

239.6° 

80.0° 

+ .5° 

270.0° 

88.3° 

256.3° 


- .2° 

125.0° 

73.2° 

186.9° 


* Only those longitude regions where 0^ 50.0° have been considered. 

** From Hakura (ref. 2) 


Table 3-4. MAXIMUM DEVIATIONS OF COMPUTED VALUES OF GEOGRAPHIC 
LONGITUDE FROM STANDARD VALUES ALONG CURVES OF 
CONSTANT GEOGRAPHIC LATITUDE * 


0 

(ao)' * 

$ 

© ** 
c 

$ ** 
c 

40.0° 

+ 

1.0° 

260.0° 

50.4° 

322.0 


- 

1.5° 

310.0° 

50.7° 

30.9 

50.0° 

+ 

1.6° 

225.0° 

52.8° 

280.5 


- 

1.1° 

340.0° 

52.3° 

65.1 

60.0° 

+ 

.6° 

240.0° 

66.3° 

290.2 


- 

1.2 

20.0 

56.4° 

99.8 

O 

o 

o 

+ 

.6° 

110.0° 

64.5° 

179.9 


- 

.7° 

25.0° 

66.1° 

109.8 

80.0° 

+ 

.2° 

310.0° 

84.4° 

90.8 

\ 

- 

.4° 

210.0° 

78.4° 

231.5 


* Only those longitude regions where Q q >_ 50.0° have been considered. 
** From Hakura (ref. 2) 
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Figure 3-1. COMPARISON OF COMPUTED (SOLID LINE) CONTOURS OF CONSTANT CORRECTED 
GEOMAGNETIC LATITUDE WITH THOSE OF HAKURA (REF. 2) (DASHED LINE) IN 
A POLAR PROJECTION OF GEOGRAPHIC COORDINATES. (COMPUTED COORDINATES 
FROM APPENDIX A) 
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In the course of determining the errors inherent in the transformations, 
it has been found that the transformations are only approximately inverses of 
one another. When geographic coordinates are converted to corrected geomagnetic 
and back again, deviations of up to 1.2 degrees latitude occur at lower latitudes. 
As with the other errors noted above, these errors decrease with increasing 
latitude. Some caution should therefore be exercised in iterative operations 
involving multiple conversions of the coordinates. 

3.3 DISCUSSION 

The corrected geomagnetic coordinates computed by Hakura (ref. 2) are based 
on the geomagnetic field for epoch 1945.0. Gustafsson (ref. 6) recently examined 
the effect of using the geomagnetic field for epoch 1965.0. He found maximum 
displacements of dipole field lines of 60 km and 70 km perpendicular to and 
along the geomagnetic meridians respectively at 0^ = 67 degrees. This corre- 
sponds to .5 degree - .6 degree geographic latitude in the context of sub- 
section 3.2 (see equation (44) and (45)). Hence, it appears that any changes 
in the coordinate system resulting from use of more recent geomagnetic field 
data will be of the same order as the errors resulting form the transformations 
presented here. It is anticipated that such changes can be adequately 
accommodated by suitable modification of the transformation parameters in 
Table 2-1. 

A basic difference between corrected geomagnetic coordinates and eccentric 
dipole coordinates, aside from the locations of the respective poles, can be 
seen by comparing the polar projection maps of Hakura (ref. 2) with those of 
Cole (ref. 7). Contours of constant eccentric dipole latitude in the latter 
are circular, whereas similar contours in the former are elliptic. This 
difference can be seen clearly in Figure 2-2. In terms of the present treat- 
ment, the transformation of ellipses to circles and vice versa is the primary 
difference between eccentric dipole coordinates and corrected geomagnetic 
coordinates in the approximation presented here. 

Gustafsson (ref. 3) has compared eccentric dipole coordinates with corrected 
geomagnetic coordinates and found the difference in latitude to be generally 
less than + 2 degrees for stations used in his study (0 >_ 61 degrees). From 

Figure 2-2, however, it appears that over certain longitude ranges deviations 
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of twice this magnitude might obtain. Thus for high latitude (0^ > 50 degrees) 
investigations in the northern hemisphere, the transformations presented here 
will likely yield better results than eccentric dipole coordinates. 

For the southern hemisphere, eccentric dipole coordinates remain a useful 
approximation. Examination of the contours of constant corrected geomagnetic 
latitude for the southern hemisphere shows that to the extent that they are 
elliptic, they are less eccentric (more nearly circular) than the corresponding 
contours for the northern hemisphere. Since, as was noted above, the major 
difference between eccentric dipole coordinates and corrected geomagnetic 
coordinates as approximated here is the ellipicity of the latitude contours, 
improvements resulting from extension of the present treatment to the southern 
hemisphere might be slight, if extant. The principal advantages of such an 
extension would then be uniformity of method for transforming coordinates in 
both hemispheres and the convenience of simple coordinate conversion in both 
directions . 

Empirical determinations in this report have been based on the 70 degree 
corrected geomagnetic latitude contour; similarly all scalings have been made 
from this contour. This CGC latitude was selected in order to minimize error 
in the vicinity of the auroral oval and because it was intermediate in the 
high latitude region. Since many investigations of high latitude phenomena 
require statistical treatment, uncertainties of 1 - 2-degrees latitude are not 
uncommon, as, for example, in the determination of the boundaries of the auroral 
oval (ref. 5). Hence, the errors determined in subsection 3.2 should not 
cause serious problems if the transformations presented here were employed in 
such studies. The simplicity and flexibility of these transformations should 
make them particularly useful in model studies, where great accuracy is not 
required, but where it is desirable to maintain the basic geometry on a large 
scale. For high precision, however, a more elaborate system such as used by 
Hakura (ref. 2) for preparing his tables or one based on even more terms in 
the harmonic expansion of the geomagnetic field may be required. 
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Section IV 
SUMMATION 


A mathematical model has been formulated for northern hemisphere contours 
of constant corrected geomagnetic latitude in a polar projection of the 
geographic coordinate system. Based on this model, analytic equations have been 
developed for transforming geographic to corrected geomagnetic coordinates 
and vice versa. These equations have been programmed for a small computer; 
programs and computed tables have been presented. 

From comparisons of computed coordinates and the coordinates on which this 
treatment was based, it has been determined that errors generally increase with 
decreasing latitude. At 50-degrees latitude maximum deviations are 2-degrees 
- 2.5-degrees latitude, improving with increasing latitude. The transformations 
are only approximately inverses of one another, so that caution must be exercised 
in iterative coordinate conversions. These transformations should be sufficiently 
accurate for use in many investigations of high latitude geophysical phenomena. 

The principal advantages of the transformations presented here are: (1) 

they are more accurate than the dipole approximations; (2) they are fast and 
easily computed on a small computer; and (3) coordinates may be converted in 
both directions between geographic and corrected geomagnetic coordinates, 
while many transformations are limited to the geographic to geomagnetic 
conversion. 
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Appendix A 

PROGRAM FOR TRANSFORMING CORRECTED GEOMAGNETIC COORDINATES 

TO GEOGRAPHIC COORDINATES 

This appendix contains the computer program MC2GC, which transforms 
corrected geomagnetic coordinates (ref. 2) to geographic coordinates. Basic 
equations for this program are developed in subsection 2.2. Characteristics 
of the program are discussed in subsection 3.1. Also presented are tables 
of computed values of geographic coordinates obtained through the input option 
for internal generation of corrected geomagnetic coordinates. 

Variables and parameters in the program are related to those in Section 
II as follows: 

GMLT = 0^ (corrected geomagnetic latitude) 

GMLG = (corrected geomagnetic longitude) 

GGLT = 0 (geographic latitude) 

GGLG = $ (geographic longitude) 

DELC = A 
RC = R 
ALPHC = a 
BETC = 6 
EPC = e 
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LOG DRIVE CART SPEC CART AVAIL PHY DRIVE 
OOOO 0001 0001 0000 

V2 M07 ACTUAL 8K CONFIG 8K 

// FOR 

•IOCS ( 1132PRINTER*CARD ) 

•ONE WORD INTEGERS 

•LIST SOURCE PROGRAM 

C MC2GC TRANSFORMS CORRECTED GEOMAGNETIC COORDINATES TO GEOGRAPHIC 

C COORDINATES BY MEANS OF AN EMPIRICAL* ANALYTIC SET OF EQUATIONS. 

C VALID ONLY FOR HIGH LATITUDES IN THE NORTHERN HEMISPHERE. 

C 

C DEFINE FUNCTIONS TO BE USED 

AMAG ( X * Y »Z ) * SQR T ( X*X + Y*Y + 2 .0*X*Y*COS ( Z ) ) 

RAD ( X ) = .0 1 74533*X 
RADI (X)=57.2957795*X 

ELPS(X.Y.Z) * SORT ((1.0 - X*X*COS(Y)*COS(Y) ) / ( 1.0 - X*X*COS ( Z ) *COS 
1 (Z ) ) ) 

C READ TRANSFORMATION PARAMETERS 

READ (2.77) DELC. RC, BETC. ALPHC. EPC 
WRITE (3.80) 

WRITE! 3.88) DELC . RC .BE TC . ALPHC . EPC 
DELC=RAD (DELC ) 

RC*RAD(RC) 

BETC*RAD (BETC) 

ALPHC=RAD( ALPHC ) 

C INPUT OPTION DESIGNATED BY SWITCH VARIABLE 10 
READ (2.78) 10 
IF (10) 11.21.11 

C OPTION TO READ IN CORRECTED GEOMAGNETIC COORDINATES (N IS THE NUMBER 

C OF LATITUDE-LONGITUDE PAIRS) 

11 READ (2.79) N 
WRITE (3.81) 

WRITE (3.82) 

DO 19 K = l.N 

READ (2.76) GMLT. GMLG 

GO TO 71 

15 WRITE (3.83) GMLT. GMLG. GGLT . GGLG 
19 CONTINUE 
CALL EXIT 

C OPTION TO GENERATE CORRECTED GEOMAGNETIC COORDINATES FOR SYNOPTIC 

C STUDIES 

21 GMLT = 45. 0 
DO 30 1=1.8 
GMLT-GMLT+5.0 
WRITE(3»81 ) 

WR I TE ( 3 » 84 ) 

WRITE(3.85) GMLT 
WR I TE ( 3 . 86 ) 

GMLG=-10 .0 
DO 29 J= 1 .36 
GMLG=GMLG+ 10.0 
GO TO 71 

25 WRITE (3.87) GMLG .GGLT .GGLG 
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29 CONTINUE 

30 CONT INUE 

call exit 

c application of transformation equations 

71 R2-RAD(90.0-GMLT) 

T 2 * R AD { GML G ) 

T1-T2-BETC 

T*T1-.13*ELPS<EPC*T1*ALPHC )*SIN(2.0*T1 ) 

R=R2*ELPS(EPC#ALPHC*T1 ) 

$ * S I N ( T ) 

C“RC/R-COS ( T ) 

T2*ATNR ( S *C ) 

I F { T 2 ) 72*73*73 

72 T2«T2+6.2831853 

73 T2»3.l4l5927+BETC-DELC-T2 
I F ( T 2* ) 74*75*75 

74 T2-T2+6. 2831853 

75 R2«AMAG(-RC*R*T) 

GGLT**90*0— RADI (R2) 

GGLG-RADI ( T 2 ) 

IF (10) 15*25*15 
C INPUT-OUTUT FORMATS 

76 FORMAT (20X*2F10.5) 

77 FORMAT ( 5 F 1 0 • 5 ) 

78 FORMAT (12) 

79 FORMAT ( 13) 

80 FORMAT ( 1 HO * 8 X * 4HDELC * 1 7X * 2 HRC , 1 7X * 4HBE TC * 1 6X • 5HALPHC * 1 6 X * 3HEPC ) 

81 FORMAT ( 1H1*26X * 7 7HTRANSF0RMAT I ON OF CORRECTED GEOMAGNETIC COORDIN 
1ATES TO GEOGRAPHIC COORDINATES) 

82 FORMAT ( 1H0*13X * 2 OHGEOMAGNE T I C LAT I TUDE t 5X *2 1HGE0MAGNET I C LONGITUD 
1E*5X* 19H GEOGRAPHIC L AT I T UDE * 6X * 2 OHGEOGRAPH I C LONGITUDE) 

83 FORMAT ( 1H * 1 3 X * 4 ( 4X * F 1 0 . 3 * 1 1 X ) ) 

84 FORMAT ( 1H0*50X • 2 2HGE0MAGNET IC LATITUDE =) 

85 FORMAT ( 1H+*73X*F10.3 ) 

86 FORMAT ( 1H0 * 2 5 X * 2 1 HGEOMAG NE T I C LONGI TUDE * 5X * 19HGE0GRAPH I C LATITUDE 
1 *5X *20HGEOGRAPHIC LONGITUDE) 

87 FORMAT ( 1 H * 2 5 X ♦ 3 ( 4X . F 1 0 . 3 * 1 1 X ) ) 

88 FORMAT ( 1H0*5(5X*F10.5*5X ) ) 

END 

FEATURES SUPPORTED 
ONE WORD INTEGERS 
IOCS 

CORE REQUIREMENTS FOR 

COMMON 0 VARIABLES 62 PROGRAM 730 

vEND OF COMPILATION 

// XEQ 


DELC 


RC 


BETC 


ALPHC 


170*00003 


EPC 
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9*50000 
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TRANSFORMATION OF CORRECTED GEOMAGNETIC COORD I NATES TO GEOGRAPHIC COORDINATES 


GEOMAGNETIC LATITUDE * 50*000 


GEOMAGNETIC LONGITUDE 

0.000 

10.000 

20.000 

30.000 
AO. 000 

50.000 

60.000 

70.000 

80.000 

90.000 

100.000 

1 10.000 

120.000 

130.000 
1 AO • 000 

150.000 

160.000 

170.000 

180.000 

190.000 

200.000 

210.000 

220.000 

230.000 
2AO.OOO 

250.000 

260.000 

270.000 

280.000 

290.000 

300.000 

310.000 

320.000 

330.000 
3 AO • 000 

350.000 


GEOGRAPHIC LATITUDE 
3 A • 2 86 
35.537 
37.376 
39.562 
A 1 • 9 1 0 
AA .310 
A6.699 
A9.021 
51.190 
53.076 
5 A • 5 28 
55 • A 1 A 
55.677 
55.365 
5A.6A8 
53.791 
53.103 

52.839 
53.103 
53.791 
5A.6A8 
55.365 
55.677 
5 5 • A 1 A 
5 A • 5 28 
53.076 
51.190 
A9.021 
A6.699 
AA. 310 
A 1 . 9 1 0 
39.562 
37.376 
35.537 
3 A . 2 86 

33.839 


GEOGRAPHIC LONGITUDE 
286. A72 
293.051 
299.881 
307.152 
315.075 
323. 83A 
333. 5A5 
3AA.23A 
355.830 
8.157 
20 • 9 3 A 
33.802 
A6.375 
58 • 3 A 8 
69.587 
80. 1A2 
90.199 

100.000 
109.800 
119.857 
130. A12 
1 A 1 . 6 5 1 
153. 62A 
166.197 
179.065 
191. 8A2 
20A. 169 
215. 76A 
226. A5A 
236.165 
2AA.92A 
252. 8A7 
260.118 
266 • 9 A 8 
273.527 

280.000 
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transformation of corrected geomagnetic coordinates to geographic coordinates 

GEOMAGNETIC LATITUDE * 55.000 

GEOMAGNETIC LONGITUDE GEOGRAPHIC LATITUDE GEOGRAPHIC LONGITUDE 


0.000 

40.070 

286. 

319 

10.000 

41.187 

292. 

736 

20.000 

42.835 

299. 

390 

30.000 

44.807 

306. 

465 

40.000 

46.944 

314. 

167 

50.000 

49.156 

322. 

680 

60.000 

51.394 

332. 

127 

70.000 

53.611 

342. 

550 

80.000 

55.731 ) 

353. 

905 

90.000 

57.631 

6. 

055 

100.000 

59.165 

18. 

767 

110.000 

60.197 

31. 

713 

120.000 

60.656 

44. 

508 

130.000 

60.571 

56. 

811 

140.000 

60.085 

68. 

435 

150.000 

59.431 

79. 

386 

160.000 

58.884 

89. 

827 

170.000 

58.672 

100. 

000 

180.000 

58.884 

110. 

172 

190.000 

59.431 

120. 

613 

200.000 

60.085 

131. 

564 

210.000 

60.571 

143. 

188 

220.000 

60.656 

155. 

492 

230.000 

60.197 

168. 

286 

240.000 

59.165 

131. 

232 

250.000 

57.631 

193. 

944 

260.000 

55.731 

206. 

094 

270.000 

53.611 

217. 

449 

280.000 

51.394 

227. 

872 

290.000 

49.156 

237. 

319 

300.000 

46.944 

245. 

832 

310.000 

44.807 

253. 

534 

320.000 

42.835 

260. 

609 

330.000 

41.187 

267. 

263 

340.000 

40.070 

273. 

680 

350.000 

39.672 

280. 

000 


\ 
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TRANSFORMATION OF CORRECTED GE O M AG N E T I C ..COORD I NAT ES TO G EOG R A PH I Q. COORD J N A T E S. 

GEOMAGNETIC LATITUDE ■ 60.000 

GEOMAGNETIC LONGITUDE GEOGRAPHIC LATITUDE GEOGRAPHIC LONGITUDE 


0,000 

45.854 


286.125 

10.000 

46.836 


292,339 

20.000 

48.291 


298.771 

30,000 

50.044 


305.598 

40,000 

51.965 


313.019 

50.000 

53.982 


321.216 

60.000 

56.059 


330.318 

70.000 

58.161 


340.382 

80.000 

60.222 


351.395 

90.000 

62.129 


3.271 

100.000 

63.742 


15.839 

1 10.000 

64.925 


28.830 

120.000 

65. 591 


41.881 

130.000 

65.747 


54.617 

140.000 

65.505 


66.774 

150.000 

65.064 


78.288 

160.000 

64.664 


89.285 

170.000 

64.504 


100.000 

180.000 

64.664 


110.714 

190.000 

65.064 


121.711 

200.000 

65.505 


133.226 

210.000 

65.747 


145.382 

220.000 

65.591 


158.118 

230.000 

64.925 


171.169 

240.000 

63.742 


184.160 

250.000 

62. 129 


196.728 

260.000 

60.222 


208.604 

270.000 

58.161 


219.617 

280.000 

56.059 


229.681 

290.000 

53.982 


238.783 

300.000 

51.965 


246.980 

310.000 

50,044 


254.401 

320.000 

48.291 


261.228 

330.000 

46.836 


267.660 

340.000 

45.854 


273.874 

350.000 

45.504 


280.000 
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transformation of corrected geomagnetic coordinates to geographic coordinates 

GEOMAGNETIC LATITUDE = 65.000 


GEOMAGNETIC longitude 

0.000 

10.000 

20.000 

30.000 

40.000 

50.000 

60.000 

70.000 

80.000 

90.000 

100.000 

110.000 

120.000 

130.000 

140.000 

150.000 

160.000 

170.000 

180.000 

190.000 

200.000 

210.000 

220.000 

230.000 

240.000 

250.000 

260.000 

270.000 

280.000 

290.000 

300.000 

310.000 

320.000 

330.000 

340.000 

350.000 


GEOGRAPHIC LATITUDE 
51.637 
52.481 
53.739 
55.268 
56.965 
58.776 
60.680 
62.651 
64.635 
66.536 
68.223 
69.560 
70.449 
70.868 
70.892 
70.682 
7 0 .440 

70.337 
70.440 
70.682 
70*892 
70.868 
70.449 
69.560 
68.223 
66.536 
64.635 
62.651 
60.680 
58.776 
56.965 
55.268 
53. 739 
52.481 
51.637 

51.337 


GEOGRAPHIC LONGITUDE 
285.872 
291.822 
297.966 
304.469 
311.522 
319.301 
327.936 
337.498 
348.004 
359.428 
11.691 
24.626 
37.938 
51.239 
64.171 
76.551 
88.422 

100.000 
111.577 
123.448 
135.828 
148.760 
162.061 
175.373 
188.308 
200.571 
211.995 
222.501 
232.063 
240.698 
248.477 
255.530 
262.033 
268.177 
274.127 

280.000 
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transformation of corrected geomagnetic coordinates TO GEOGRAPHIC COORDINATES 

GEOMAGNETIC LATITUDE = 70.000 


geomagnetic longitude 

GEOGRAPHIC LATITUDE 

GEOGRAPHIC LONGITUDE 

0,000 

57.419 

265.531 

10.000 

58.123 

291.12.3 

20.000 

59.178 

296.876 

30.000 

60.475 

302 .94 2 

AO. 000 

61.935 

309. 494 

50.000 

63.523 

316.696 

60.000 

65.230 

324.671 

70.000 

67.043 

333.495 

80.000 

68.923 

343.209 

90.000 

70.789 

353.849 

100.000 

72.531 

5.450 

110.000 

74.022 

18.019 

120.000 

75.154 

31.446 

130.000 

75.874 

45.434 

140.000 

76.209 

59.545 

150.000 

76.268 

73.395 

160.000 

76.208 

86.838 

170.000 

76.169 

100.000 

180.000 

76.208 

113.161 

190.000 

76.268 

126.604 

200.000 

76.209 

140.454 

210.000 

75.874 

154.565 

220.000 

75.154 

168.553 

230.000 

74.022 

181.980 

240.000 

72.531 

194.549 

250.000 

70.789 

206.150 

260.000 

68.923 

216.790 

270.000 

67.043 

226.504 

280.000 

65.230 

235.328 

290.000 

63.523 

243.303 

300.000 

61.935 

250.505 

310.000 

60.475 

257.057 

320.000 

59.178 

263.123 

330.000 

58.123 

268.876 

340.000 

57.419 

274.469 

350.000 

57.169 

280.000 
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transformation of corrected geomagnetic coordinates to geographic coordinates 

GEOMAGNETIC LATITUDE = 75.000 

GEOMAGNETIC LONGITUDE GEOGRAPHIC LATITUDE GEOGRAPHIC LONGITUDE 


0.000 

63.199 

285.041 

10.000 

63.768 

290.122 

20.000 

64.602 

295.321 

30.000 

66.661 

300.763 

40.000 

66.854 

306.598 

50.000 

68.192 

312.965 

60.000 

69.665 

319.961 

70.000 

71.272 

327.644 

BO. 000 

72.989 

336.050 

90.000 

74*759 

345.240 

100.000 

76.500 

355.349 

1 10.000 

78.112 

6.608 

120.000 

79.499 

19.294 

130.000 

80.582 

33.577 

140.000 

81.326 

49.312 

150.000 

81*755 

65.982 

160.000 

81.950 

82.986 

170.000 

82.002 

100.000 

180.000 

81.950 

117.013 

190.000 

81.755 

134.017 

200.000 

81.326 

150.687 

210.000 

80.582 

166.422 

220.000 

79.499 

180.706 

230.000 

78.112 

193.391 

240.000 

76. $00 

204.650 

250.000 

74.759 

214.759 

260.000 

72.989 

223.949 

270.000 

71.272 

232.355 

280.000 

69.665 

240.039 

290.000 

68.192 

247.034 

300.000 

66.854 

253.401 

310.000 

65.651 

259.236 

320. 000 

64.602 

264.678 

330.000 

63.758 

269.877 

340.000 

63. 199 

274.958 

350.000 

63.002 

280.000 
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transformation of corrected geomagnetic coordinates to geographic coordinates 


GEOMAGNETIC LONGITUDE 
0.000 
10.000 
20.000 

30.000 

40.000 

50.000 

60.000 

70.000 

80.000 

90.000 

100.000 

110.000 

120.000 

130.000 

140.000 

150.000 

160.000 

170.000 

180.000 

190.000 

200.000 

2 10.000 

220.000 

230.000 

240.000 

250.000 

260.000 

270.000 

230.000 

290.000 

300.000 

3 10.000 

320.000 

330.000 

340.000 

350.000 


GEOMAGNET I 

C LATITUDE 

GEOGRAPHI 

C LATITUDE 

68 

.976 

69 

.380 

69 

.996 

70 

.773 

71 

.683 

72 

.718 

73 

.887 

75 

.195 

76 

.630 

78 

.159 

79 

.731 

81 

.293 

82 

.797 

84 

.208 

85 

.498 

86 

.626 

87 

.485 

87 

.834 

87 

.485 

86 

.626 

85 

.498 

84 

.208 

82 

.797 

81 

.293 

79 

.731 

78 

.159 

76 

.630 

75 

.195 

73 

.887 

72 

.718 

71 

.683 

70 

.773 

69 

.996 

69 

.380 

68 

.976 

68 

.834 


80.000 


GEOGRAPHIC 

LONGITUDE 

284 

• 283 

288 

.575 

292 

.924 

297 

.415 

302 

.159 

307 

.242 

312 

.710 

318 

.549 

324 

.719 

331 

.185 

337 

.989 

345 

.318 

353 

.579 

3 

.503 

16 

.336 

34 

. 308 

61 

.367 

100 

.000 

138 

.632 

165 

.691 

183 

.663 

196 

.496 

206 

.420 

214 

.681 

222 

.010 

228 

.814 

235 

.280 

241 

.450 

247 

.289 

252 

.757 

257 

.840 

262 

.584 

267 

.076 

271 

.424 

275 

.716 

280 

.000 
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TRANSFORMATION OF CORRECTED GEOMAGNETIC COORDINATES TO GEOGRAPHIC COORDINATES 

GEOMAGNETIC LATITUDE = 85.000 


GEOMAGNETIC LONGITUDE 

0.000 

10.000 

20.000 

30.000 

40.000 

50.000 

60.000 

70.000 

80.000 

90.000 

100.000 

110.000 

120.000 

130.000 

140.000 

150.000 

160.000 

170.000 

180.000 

190.000 

200.000 

210.000 

220.000 

230.000 

240.000 

250.000 

260.000 

270.000 

280.000 

290.000 

300.000 
3 10.000 

320.000 

330.000 

340.000 

350.000 


GEOGRAPHIC LATITUDE 
74.747 
74.976 
75.330 
75.783 
76.324 
76.953 

77.676 
78.495 

79.399 
80.363 
81.352 
82.333 
83.276 
84. 163 
84.973 
85.662 
86.151 
86.332 
86.151 
85.662 
84.973 
84.163 
83.276 
82. 3'33 
81.352 
80.363 

79.399 
78.495 

77.676 
76.953 
76.324 
75.783 
75.330 
74.976 
74.747 
74.667 


GEOGRAPHIC LONGITUDE 
282.950 
285.873 
288.771 

291.676 
294.631 
297.648 
300.687 
303.639 
306.341 
308.600 

310.228 
311.058 
310.916 
309.537 
306.424 
300.760 

291.770 
• 280.000 

268.229 
259.239 
253.576 
250.462 
249.083 
248.941 

249.771 

251.399 
253.659 
256.360 
259.312 
262.351 
265.368 
268.323 

•• 271.228 

274.126 
277.049 

280.000 
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Appendix B 

PROGRAM FOR TRANSFORMING GEOGRAPHIC COORDINATES 
TO CORRECTED GEOMAGNETIC COORDINATES 

This appendix contains the computer program GC2MC, which transforms 
geographic coordinates to corrected geomagnetic coordinates (ref. 2). Basic 
equations for this program are developed in subsection 2.3. Characteristics 
of the program are discussed in subsection 3.1. Also presented are tables of 
computed values of corrected geomagnetic coordinates obtained through the input 
option for internal generation of geographic coordinates. 

Variables and parameters in the program are related to those in Section 
II as follows: 

GMLT = 6^ (corrected geomagnetic latitude) 

GMLG - (corrected geomagnetic longitude) 

GGLT = 0 (geographic latitude) 

GGLG = $ (geographic longitude) 

DELC = A 
RC = R 
ALPHC = a 
BETC = 6 
EPC = e 
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LOG DRIVE CART SPEC CART AVAIL PHY DRIVE 
OOOO 0001 0001 0000 

V? M 0 7 ACTUAL 8 K “CONF I G “ 8K 

// FOR 

*IOCS( 1132PR INTER .CARD) 

♦ONE WORD INTEGERS 

♦LIST SOURCE PROGRAM _ _ 

C GC2MC TRANSFORMS GEOGRAPHIC COORD I NATES ' TO CORRECTED GEOMAGNETIC 
C COORDINATES BY MEANS OF AN EMPIRICAL* ANALYTIC SET OF EQUATIONS. VALID 
C ONLY FOR HIGH LATITUDES IN THE NORTHERN HEMISPHERE. 

C ONLY FOR HIGH LATITUDES IN THE NORTHERN HEMISPHERE. 

C DEFINE FUNCTIONS TO BE USED 

AMAG(X,Y*Z> = SORT ( X*X + Y*Y + 2 .0*X^Y^COS ( ZJ ) _ 

RAD ( X ) = . 0 1 745 33*X 
RADI ( X )=57.2957795*X 

ELPS( X*Y *Z ) = SORT ((1.0 - X*X*COS( Y) *COS ( Y ) ) / ( 1 . 0 - X*X*COS < Z > *COS 
1 ( Z ) ) ) 

C READ TRANSFORMATION PARAMETERS 

READ (2,77) DELC. RC ♦ BETC. ALPHC • EPC 

WRITE! 3.80) 

WRI TE< 3.8 8 ) DELC. RC. BETC. ALPHC. EPC 
OELC=RAD ( DELC ) 

RC=RAD ( RC ) 

BETC=RAD (BETC ) 

ALPHC=RAD( ALPHC ) 

C INPUT OPTION DESIGNATED BY SWITCH VARIABLE 10 

READ (2,78) 10 

IF (10) 11.21.11 

C OPTION TO READ IN GEOGRAPHIC COORDINATES ( N IS NUMBER OF LATITUDE- 
C LONGITUDE PAIRS) 

11 READ (2.79) N 
WRITE (3.81) 

WRITE (3.82) 

DO 19 K = l.N 

READ (2.77) GGLT . GGLG 

GO TO 71 

15 WRITE (3.83) GMLT . GMLG . GGLT. GGLG 
19 CONTINUE 
CALL EXIT 

C OPTION TO GENFRATE GEOGRAPHIC COORDINATES FOR SYNOPTIC STUDIES 
21 GGLT = 45.0 
DO 30 1=1.8 
GGLT-GGLT + 5.0 

WRITE (3. 81) - - - 

WRI TE ( 3.84 ) 

WRITE (3.85) GGLT 
WRITE (3.86) 

GGLG = — 10 .0 
DO 29 J= 1 » 36 
GGL G = GGLG+ 10 .0 
GO TO 71 

25 WR I TE ( 3 » 87 ) GGLG .GMLT .GMLG 
29 CONTINUE 
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30 CONTINUE 
call EXIT 

C APPLICATION OF TRANSFORMATION EQUATIONS 
71 R1=RAD(90.0 - GGLT) 

T 1 =R AD ( GGLG) _ __ 

ti=ti-betc+delc ~ 

S=S IN ( T1 J 
C=RC/R 1 + COS ( T 1 ) 

T = ATNR ( S *C ) 

T3 = T + • 1 3 *E L PS ( EPC * T * AL PHC ) *$ I N ( 2 • 0 *T ) 

I F ( T 3 ) 73*75*75 

73 T3 = T3 + 6.2831S53 “ ™ ' 

7 5 R = AMAG(RC *R1*T1 ) *ELPS ( EPC * T . ALPHC J 
GMLTs90.0-RADI ( R > 

GMLGaRADI ( T3+BETC ) 

IF ( G M LG-3 60 * 0 ) 62*62*61 

61 GMLG=GMLG-360*0 

62 IF (10) 15*25*15 

C INPUT-OUTUT FORMATS 

77 FORMAT ( 5F10# 5 ) 

78 FORM AT ( I 2 ) 

79 FORMAT (13) 

90 FORMAT ( 1 HO * 9 X » 4HDELC * 1 7X » 2HRC * 1 7X >4HBETC * 16X ♦ 5HALPH C > 16X * 3HEPC ) 

91 FORMAT ( 1 H 1 * 26X * 77Hf RANSFORMAT TON OF GEOGRAPHIC COORDl'NATES TO COR 
1RECTED GEOMAGNETIC COORDINATES) 

82 FORMAT ( 1H0. 13X . 19HGE0GRAPHIC L AT I T UDE * 6 X * 20HGE0GRAPH I C LONGITUDE* 
15X*20HGEOMAGNETIC LAT I TUD E . 5X . 2 1 HGEOMAGN E T I C N.ONGITUDE) 

83 F0RMAT(1H * 13X .4 ( 4X.F10.3 * UX ) ) 

84 FORMAT ( 1 HO . 50X * 2 i HGEOGRAPH I C LATITUDE * ) 

95 FORMAT (lH+* 73XiFlb.3'j “ ‘ ' ' 

. 86 FORMAT ( 1H0.25X * 2 OHGEOGRAPH I C LONG I TUDE * 5X *20HGEOMAGNET I C LATITUDE 
1 * 5X * 2 lHGEOMAGNET I C LONGITUDE) 

87 FORMAT (1H *25X*3(4X*F10.3*11X) ) 

89 FORMAT (1H0*5(5X*F10.5*5X) ) 

END 


FEATURES SUPPORTED 
ONE w’ORD INTEGERS 
IOCS 

CORE REQUIREMENTS FOR 

COMMON 0 VARI ABLES "6 2 ' ^PROGRAM 112 

end of compilation 

// XEO 


DELC 

70.00001 


RC 

9.50000 


BETC 

170.00003 


ALPHC 

55.00000 



EPC 

0.59130 
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T RANS FORMATION .0 F. GEOGRAPHIC C OORDINATES TO CQRRECTFD GFOMAGNF TIC COORDINATES 



GEOGRAPHIC LATITUDE = 

50.000 

geographic longitude 

GEOMAGNETIC LATITUDE 

geomagnetic longitude 

0*000 

48.068 

82.687 

10.000 

46 • 7 06 

89.833 

20.000 

45.816 

97.178 

30.000 

45.375 

104.936 

40.000 

45.328 

113.197 

50.000 

45.591 

121.950 

60.000 

46.054 

131.122 

70.000 

46.596 

140.609 

80.000 

47.093 

150.306 

90.000 

47.441 

160.126 

100.000 

47.565 

169.999 

110.000 

47.441 

179.873 

120.000 

47.093 

189.693 

130.000 

46.596 

199.390 

140.000 

46.054 

208.877 

150.000 

45.591 

218.048 

160. 000 

45.328 

226.802 

170.000 

45.375 

235.062 

180.000 

45.816 

242.821 

190.000 

4.6. 706 

250.166 

200.000 

48.068 

2 57.312' 

210.000 

49.883 

264.610 

220.000 

52*092 

272.531 

230.000 

54.586 

281.591 

240.000 

57.205 

292.229 

250.000 

59.726 

304.658 

260.000 

61.870 

318.750 

270.000 

63.331 

334.064 

280.000 

63.853 

350.000 

290.000 

63.331 

5.935 

300.000 

61.870 

21.249 

310.000 

59.726 

35.341 

320.000 

57.205 

47.770 

330.000 

54.586 

58.408 

340.000 

52.092 

67.468 

350.000 

49.883 

75.389 
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tra nsformation of geographic coordinates to corrected geomagnetic coordinates 



GEOGRAPHIC LATITUDE » 

55.000 

geographic longitude 

GEOMAGNETIC LATITUDE 

GEOMAGNETIC LONGITUDE 

0*000 

53.242 

84.115 

10.000 

51.913 

91.230 

20.000 

51.002 

98.556 

30.000 

50.489 

106.277 

AO. 000 

50.327 

11A.A62 

50.000 

50.443 

123.090 

60.000 

50.745 

132.089 

70.000 

51.131 

1A1.365 

80.000 

51.497 

150. 82A 

90.000 

51.758 

160.388 

100.000 

51.851 

169.999 

110.000 

51.758 

179.610 

120.000 

51.497 

189.175 

130.000 

51.131 

198. 63A 

140.000 

50.745 

207.910 

150.000 

50.443 

216.909 

160.000 

50.327 

225.537 

170.000 

50.489 

233.722 

180.000 

51.002 

241.443 

190.000 

51.913 

248.769 

200.000 

53.2A2 

255.884 

210.000 

5A.97A 

263.114 

220.000 

57.057 

270.922 

230.000 

59.398 

279.857 

2AO.OOO 

61.854 

290.439 

250.000 

64.225 

302.983 

260.000 

66.252 

317. A35 

270.000 

67.641 

333. 3A0 

280.000 

68.139 

350.000 

290.000 

67.641 

6.659 

300.000 

66.252 

22.564 

310.000 

64.225 

37.016 

320.000 

61.854 

A9.560 

330.000 

59.398 

60.1A2 

3 AO . 000 

57.057 

69.077 

350.000 

5A.97A 

76.885 
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TRANSFORMATION. QE. JjEOGlRAEHIC ..COOADINATlS TO. CORRECT LD. GEOMAGNETIC CQQRD1.NAT.ES 

GEOGRAPHIC LATITUDE = 60.000 

GEOGRAPHIC LONGITUDE GEOMAGNETIC LATITUDE GEOMAGNETIC LONGITUDE 


0.000 

58.379 

86 

.018 

10.000 

57.086 

93 

.086 

20.000 

56.157 

100 

.372 

30.000 

55.576 

108 

.027 

AO. 000 

55.303 

116 

.095 

50.000 

55.277 

124 

.547 

60.000 

55.423 

133 

.317 

70.000 

55.657 

142 

.319 

30.000 

55.897 

151 

.47 5 

90.000 

56.073 

160 

.718 

100.000 

56.138 

169 

.999 

110.000 

56.073 

179 

.281 

120.000 

55.897 

188 

. 52A 

130.000 

55.657 

197 

.680 

1A0.000 

55.423 

206 

.682 

150.000 

55.277 

215 

• 451 

160.000 

55.303 

223 

.904 

170.000 

55.576 

231 

*972 

180.000 

56.157 

239 

.626 

190.000 

' 57.086 

246 

.913 

200.000 

. 58.379 

253 

.981 

210.000 

60.026 

261 

.118 

220.000 

61.982 ■ 

268 

.768 

230.000 

64.171 

277 

.508 

240.000 

66.467 

287 

.961 

250.000 

68.696 

300 

.604 

260.000 

70.6 17 

315 

. 5 2 2 

270.000 

7.1 .946 

332 

.269 

280.000 

72.426 

, 3 50 

.000 

290.000 

71.946 

* 7 

.730 

300.000 

70.617 

24 

.477 

310.000 

68.696 

39 

.396 

320.000 

66.467 

52 

.038 

330.000 

64.171 

62 

.491 

340.000 

61.982 

71 

.231 

350.000 

60.025 

78 

.881 
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TRANS FORMAT I UN OF GEOGRAPHIC COORDINATES TO CORRECTED GEOMAGNETIC COORDINATES 

GEOGRAPHIC LATITUDE = 6S.000 


GEOGRAPHIC LONGITUDE 

GEOMAGNETIC LATITUDE 

GEOMAGNETIC LONGITUDE 

0.000 

63.457 

88.679 

10.000 

62.205 

95.665 

20.000 

61.264 

102.869 

30.000 

60.622 

110.400 

40.000 

60.244 

118.279 

50.000 

60.084 

126.474 

60.000 

60.082 

134.924 

70.000 

60.172 

143.560 

30.000 

60.292 

152.317 

90.000 

60.388 

161.143 

100.000 

60.424 

169.999 

110.000 

60.388 

178.856 

O 

o 

o 

• 

c 

f\. 

60.292 

187.682 

1 30 • 000 

60.172 

196.439 

140*000 

60.082 

205.075 

150. 001 

60.084 

213.525 

160. 00C 

60.244 

221.720 

170.000 

60.622 

229.599 

180.000 

61.264 

237.130 

190.000 

62.205 

244.334 

200.000 

63.457 

2 5 1 . 3 2 0 

210.000 

65*014 

258 .324 

220.000 

66.842 

265. 743 

230.000 

68.877 

274. 164 

240.000 

71.019 

284.330 

250.000 

73.116 

296.976 

260.000 

74.950 

312.488 

270. 000 

76.241 

330.524 

280.000 

76.712 

350.000 

290.000 

76.241 

9.475 

300.000 

74.950 

27.511 

310.000 

73.1 16 

43.023 

320.000 

71.019 

55.669 

330.000 

68.877 

65.835 

340.000 

66.842 

74.256 

350.000 

65.014 

81.675 
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M-765 


NORTHROP 

HUNTSVILLE " 


J-R ANS FORMATION OF GEOGRAPHIC COORDINATES TO CORRECTED GEOMAGNETIC COORDINATES 



GEOGRAPHIC LATITUDE - 

70.000 

GEOGRAPHIC LONGITUDE 

GEOMAGNETIC LATITUDE 

GEOMAGNETIC LONGITUDE 

0.000 

68.434 

92.650 

10.000 

67.233 

99.471 

20.000 

66.292 

106.492 

30.000 

65.600 

113.778 

40.000 

65.131 

121.335 

50.000 

64.849 

129.129 

60.000 

64.711 

137.114 

70.000 

64.669 

145.236 

80.000 

64.678 

153.448 

90 • 000 

64.700 

161.712 

100.000 

64.710 

169.999 

110.000 

64.700 

178.287 

120.000 

64.678 

186.551 

130.000 

64.669 

194.763 

140.000 

64.711 

202.885 

150.000 

64.849 

210.870 

160.000 

65.131 

218.664 

170.000 

65.600 

226.221 

180.000 

66.292 

233.507 

190.000 

67.233 

240.528 

200.000 

• 68.434 

247.349 

210.000 

69.892 

254.134 

220.00 0" 

71.583 

261.196 

230.000 

73.46-1 

269.071 

240.000 

75.451 

278.589 

250.000 

77.432 

290.870 

260.000 

79.213 

306.996 

270.000 

80.511 

327.185 

290.000 

•80.998 

350.000 

290.000 

80.511 

12.814 _ _ 

300.000 

79.213 

' 33.001 

310.000 

77.432 

49.129 

320.000 

75.451 

61.410 

330.000 

73.461 

70.928 

340.000 

71.583 

78.803 

350.000 

69.892 

85.865 
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Momsfiiop 

HUNTSVILLE " 


TRANS format ion of geographic coordinates to corrected geomagnetic coordinates 


GEOGRAPHIC LONGITUDE 
0*000 
10*000 

GEOGRAPHIC LATITUDE = 

GEOMAGNETIC LATITUDE 

"'737221 

72.095 

75.000 

GEOMAGNETIC LONGITUDE 
99.145 
105.557 

20.000 

71.178 

112 

► 1 34 

30.000 

70.460 

118 

>907 

40.000 

69.923 

125 

>871 

50.000 

69.543 

133 

>000 

60.000 

69.291 

140 

>263 

70.000 

69.136 

147 

>623 

90.000 

69.050 

155 

.049 

90.000 

69.009 

162 

.515 

100.000 

68.997 

169 

.999 

110.000 

69.009 

177 

.484 

120.000 

69.050 

184 

.950 

130.000 

• 69.136 

192 

.376 

140.000 

69.291 

199 

.736 

150.000 

69.543 

206 

► 999 

160.000 

69.923 

214 

► 128 

170.000 

70.460 

221 

► 092 

180.000 

71.178 

227 

.865 

190.000 

72.095 

234 

.442 

"200.000 

73.221 

240 

► 8 54 

210.000 

74.555 

247 

► 185 

220.000 

76.085 

253 

► 619 

230.000 

77.704 

260 

>511 

240.000 

79.609 

268 

>552 

250.000 

81.484 

279 

>090 

260.000 

83.272 

294 

>620 

270.000 

84.699 

318 

>375 

280.000 

85.285 

350 

>000 

290.000 

84.699 

21 

624 

300.000 

83.272 

45 

>379 

310.000 

81.484 

60 

>909 

" 320.000 

79.609 

71 

>447 

330.000 

77.784 

79 

488 

340.000 

76.085 

86 

3 80 

350.000 

74.555 

92 

814 
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NORTHROP 

HUNTSVILLE " 


-.-TRANSFORMATION OF QE 

ographic coordinates to corrected geomagnetic coordinates 


GEOGRAPHIC LATITUDE « 

80.000 

geographic LONGITUDE 

GEOMAGNETIC LATITUDE 

geomagnetic longitude 

0.000 

77.595 

111.085 

10.00C 

76.608 

116.309 

20.000 

75.774 

121.736 

30.000 

75.085 

127.355 

40.000 

74.531 

133.145 

50.000 

74.101 

139.080 

60.000 

73.779 

145.133 

70.000 

73.549 

151.276 

90.000 

73.397 

157.483 

90.000 

73.311 

163.731 

100.000 

73.283 

169.999 

110.000 

73.311 

176.268 

120.000 

73.397 

182.516 

130.000 

73.549 

188.723 

140.000 

73.779 

194.866 

150.000 

74.101 

200.919 

160.000 

74.531 

206.854 

170.000 

75.085 

212.644 

190.000 

75.774 

218.263 

190.000 

76.608 

223.690 

200.000 

77.595 

228.914 

210.000 

78.733 

233.936 

220.000 

80.018 

238.776 

220.000 

81.441 

243.484 

240.000 

82.987 

248. 162 

250.000 

84.636 

253.040 

260.000 

86.363 

258.750 

270.000 

88.130 

268.303 

290.000 

89.571 

350.000 

290.000 

88.130 

71.696 

300.000 

86.363 

81.249 

310.000 

84.636 

86.959 

320.000 

82.987 

91.837 

330.000 

81.441 

96.516 

340.000 

80.018 

101.223 

350.000 

78.733 

106.063 
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NORTHROP M-765 

HUNTSVILLE — — — — — ^ "" ~ 


transformation of geographic co ord inat es to correcte d geoma gnetic coord inates 

GEOGRAPHIC LATITUDE «. . . 85.000 

GEOGRAPHIC LONGITUDE GEOMAGNETIC LATITUDE GEOMAGNETIC LONGITUDE 


0.000 

80.936 

134.449 

10.000 

80.J75 

136.679 

20.000 

79.684 

139.408 

30.000 

79.169 

142.528 

40.000 

78.730 

145.956 

50.000 

78.365 

149.629 

60.000 

78.072 

153.491 

70.000 

77.849 

157.499 

80.000 

77.693 

161.610 

90.000 

77.600 

165.788 

100.000 

77.569 

169.999 

110.000 

77.600 

174.211 

120.000 

77.693 

178.389 

130.000 

77.849 

182.500 

140.000 

78.072 

186.507 

150. 00C 

78.365 

190.370 

160.000 

78.730 

194.043 

170.000 

79.169 

197.471 

180.000 

79.684 

200.591 

190.000 

80.275 

203.320 

20 0.000 

80.936 ^ 

205.550 

210.000 

81.660 

207.135 

220.000 

82.437 

207.864 

230.000 

83.246 , 

207.428 

240.000 

84.058 

205.372 

250.000 

84.832 

201.060 

260.000 

85.500 

""**193.760 

270.000 

85.970 

s 183.104 

280.000 

86.142 

169.999 

290.000 

85.970 

156.895 

300.000 

85.500 

146.239 

310.000 

84.832 

138.938 

320.000 

84.058 

'134.627 " 

330.000 

83.246 

132.571 

340.000 

82.437 

132.135 

350.000 

81.660 

132.864 
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